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Abstract 

Multiplication of spin qubits arises at double resonance in a bichromatic field when the frequency 
of the radio-frequency (rf) field is close to that of the Rabi oscillation in the microwave field, pro- 
vided its frequency equals the Larmor frequency of the initial qubit. We show that the operational 
multiphoton transitions of dressed qubits can be selected by the choice of both the rotating frame 
and the rf phase. In order to enhance the precision of dressed qubit operations in the strong-field 
regime, the counter-rotating component of the rf field is taken into account. 

PACS: 03.67.Pp, 33.40.+f, 33.35. -t-r 

Theoretical models of quantum computations assume the existence of an ideal two-level quan- 
tum system (qubit) and the possibility of an exact description of the qubit's interaction with 
external electromagnetic fields [1] . It is known that the resonant interaction between electromag- 
netic radiation and qubit induces Rabi oscillations, which are the basis for quantum operations. 
The Rabi frequency ur is defined by the amplitude of the electromagnetic field and usually 
is much smaller than the energy difference Uq (in frequency units) between the qubit's states. 
The "dressing" of qubit by the electromagnetic field splits each level into two giving rise to two 
new qubits with energy difference ur. The spectrum of the multilevel "qubit -|- field" system 
consists of three lines at the frequencies ojq and ojq ± ojr (the Mollow triplet [2]). The second 
low-frequency electromagnetic field with the frequency close to the Rabi frequency ujr could in- 
duce an additional Rabi oscillation on dressed states of new qubits. These qubits are attracting 
interest because their coherence time is longer than that of the initial qubit [3-5]. The results 
of studies of qubits dressed by bichromatic radiation formed by fields with strongly different 
frequencies are important for a wide range of physical objects, including, among others, nuclear 



and electron spins, double-well quantum dots, flux and charge qubits in superconducting sys- 
tems. In NMR [6, 7], EPR [5, 8, 9] and optical resonance [10] such investigations are used in the 
development of line-narrowing methods. 

In this letter, we describe the multiplication of spin qubits at double resonance in a bichro- 
matic field with strongly different frequencies. We then show that the operational multiphoton 
transitions of dressed qubits can be selected by the choice of both the rotating frame and the 
phase of the low-frequency field. Two important examples of such transitions in the rotating and 
doubly rotating frames are presented. 

Let an electron spin qubit be in three fields: a microwave (mw) one directed along the x 
axis of the laboratory frame, a radio-frequency (rf) one directed along the z axis, and a static 
magnetic one also directed along the z axis. The Hamiltonian of the qubit in these fields can be 
written as follows: 

H = Ho + H^ (t) + (t) . (1) 

Here Hq = uqS^ is the Hamiltonian of the Zeeman energy of a spin in the static magnetic field Bq, 
where Uq = '^Bq, and 7 is the electron gyromagnetic ratio. Moreover, (t) = 2uJi cos {uot + 0) 
and H\\ (t) = 2uj2 cos [urft + ip) are the Hamiltonians of the spin interaction with linearly po- 
larized mw and rf fields, respectively. Bi and B2, OJ and cu^/, and ip and ifj denote the respective 
amplitudes, frequencies, and phases of the mw and rf fields. Finally, ui = ■jBi and 0J2 = 7-B2 
stand for the Rabi frequencies, whereas s^'^'^ are the components of the spin operator. 

The evolution of the system with the Hamiltonian [T] is described by the Liouville equation 
for the density matrix p: 



(2) 



(we set the Planck constant h = 1). We perform the transformation (p — pi = UipUi, 
Ui = e"*'^*'*^) to the singly rotating frame, which rotates with frequency u around the z axis 
of the laboratory frame. In this frame, Eq. [2] turns into: 

^^ = [HuPl], (3) 

where Hi = UiHUi = As^ + {ui/2){s'^ + s') + 2uj2 cos{urft + iIj)s^ and A = uq — u. The mw 
phase = and the counter-rotating component of the mw field is neglected. We also assume 
that the exact resonance condition is fulfilled A = 0, and that Ui, Urf ^ 002- Upon rotation 
of the frame around the y axis by the angle of 7r/2 (pi — > p2 = U2P1U2, U2 = e"*'^'^''/^, where 
= (s+ — s^)/2i), we obtain: 

t^ = [H2,P2], (4) 

where H2 = U2H1U2 = Wis^ — 002 cos{urft + 4'){s~^ + s^). 

Now, we pass to the interaction representation by choosing the frame rotating with frequency 
Ui around the z axis (p2 Ps = P2U3, f/3 = e"*'^^*'*''). In this frame we have: 

^^ = [^3,P3], (5) 

where 

= U^H2Ui = -{uj2/2)s+ (e*'5*e-*'^ + e*(2-^/+'^)*e*'^) - h.c, 6 = ooi - oorf, and \6\ < Wi, cOrf 
in our case. Rapidly oscillating (e^*^'^'^^*) terms in the Hamiltonian can be eliminated by 
the Krylov-Bogoliubov-Mitropolsky method [5, 11, 12]. Averaging over the period 2'ir/ujrf, we 
obtain the following effective Hamiltonian up to the second order in uj2/ujrf'- 



H, ^ H^ff = + Z/^) . (6) 



In the above equation we have put: 

i^g^ = {H,{t)) = -{oo,/2){sU^'^e-^^ + h.c), = ^^{[fdrmr)-mr))),H,m = 

Abss^. The symbol < ... > denotes time averaging: {A{t)) = ^ A{t)dt, where T = 27T/uJrf 
and Ab5 ~ a;|/4a;r/ is the Bloch-Siegert-hke frequency shift. 

After the canonical transformation Ps ^ Pi = p^iU^^, f/4 = ^-^i^^^-i')^'' ^ the equation 



^^ = [^e//,P3] (7) 

is transformed into 



^^ = [^4,P4], (8) 

where = U^HeffUi = {6 + Abs)s' - (w2/2)(s+ + s"). 

The diagonalization of the Hamiltonian H4 by means of the rotation operator U5 = e~*^*^ 
(P5 ^ P5 = U^pdJ^, = U^HdJ^) yields: 

^^ = [^5,P5]. (9) 

Here = es^, e = a/ {ooi — ujrf + AbsY + ^2 the frequency of the Rabi oscillations be- 
tween the spin states dressed simultaneously by the mw and rf field while sin,^ = —1^2/ s, 

COS^ = (Wi - Urf + Abs)/£- 

By using Eqs. [2] - [9l the density matrix in the laboratory frame (LF) can be written as: 

-iH5t„{n\JH5tTT+ 



p{t) = t/if/2f/3t/4f/5e-*^^V5(0)e*^^*f/5+t/+t/3+t/2+f/+, (10) 



where 



P5(0) = f/5+f/+(0)f/3+(0)[/2+f/+(0)p(0)f/5[/4(0)f/3(0)[/2t/i(0), (11) 

and f/i(0) = 1, t/3(0) = 1, ?74(0) = e-'^'\ 
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Figure 1: Energy-level diagram of a qubit and transitions created by a bicliromatic field at double 
resonance {uq = u, Ui = Urf). 

Initially, the qubit is in the ground state and p(0) = 1/2 — s^. The absorption signal 
v(t) = Sp{p(t){s^ — s^)}/2i in the laboratory frame can be derived from Eqs. [TOlandfTT] 

v{t) = {{l\p{t)\2) - {2\p{t)\l)) /2z = 

= (1/2) sin,^ cos,^ cosip smut + (1/16){4 sin,^ simp[cos{uj + e)t — cos{u — e)t] + 
— 4sin^ cos^ COS'?/' [sin — e)t + sm{u + e)t] + 2sin^^ sin 2-?/; [cos (a; + 0Jrf)t + cos^u — Urf)t] + 
+ ((cos,^ — 1)^ + (cos^,^ — l) cos2?/') X [sin(a; + Urf — e)t — sin(a; — a;,./ + s)t\ + 
+ (cos^^ — l) sin 2?/) X [cos(a; — Urf + e)t + cos(a; + Urf — £)t\ + 



((cos ^ + 1)^ + (cos^ ^ ~~ l) cos 2il)) X [sin(ci; + Urf + e)t — sin(c<; — Urf — + 
+ (cos^,^ — l) sin 2'?/' X [cos(ci; — cUr/ — s)t + cos(a; + w^/ + 



(12) 



The resonant interaction between the mw field and the qubit creates its dressed states and 
two new qubits with energy splitting equal to the Rabi frequency Wi, as shown in Fig.l. The rf 



field with the frequency cu^/, which is close to the Rabi frequency uj\ of the new qubits, "dresses" 



these qubits, giving rise to four qubits with the energy sphtting e. Allowed transitions between 
states of these qubits afford nine spectral lines observed in the laboratory frame. 

Figs. 2 and 3 show the time evolution of absorption signals and their Fourier spectra of 
dressed qubits under conditions typical for EPR. 
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Figure 2: Time evolution of the absorption signals Figure 3: Fourier spectra of the absorption signals in 

in the laboratory (solid line), singly rotating (dot the laboratory frame shown in Fig. 2 by a solid line. 

line) and doubly rotating (dash line) frames. The 

signals were obtained for the following parameters of 

the bichromatic field: to = ujq, wi = w^/ = SttI.O 

MHz, U2 = 27r0.24 MHz, using the exponential decay 

function with T = 16 fis. 



For the rf phase ip = 0, three triplets with the intensive central lines at w — cJi, cj and u + uJi 
are formed (Fig. 3a). The less intensive sidebands have the frequencies ±e relative to each of 
the central lines. For the rf phase ifj = 7r/2, the central lines in the triplets vanish (Fig. 3b), 
each triplet turning into an doublet. When the rf phase is random, averaging over a sufficiently 
large number of experiments (at the uniform distribution of the phase in the interval from to 
2n) leads to a complete removal of the central triplet. The differences of line's intensities in the 
two residual triplets {u — Ui, u — Ui ± e and u + Ui, u + Ui ± e) become smaller (Fig. 3c). 

There is the possibility of selecting the observed transitions of four qubits by employing the 
rotating frame. In the singly rotating frame (SRF), the absorption signal described by the density 
matrix pi{t) (Eq. (4)) can by written as 

vi (t) = (1/8) [2 sin^ ^ (sin w^/t + sin {urft + 2iIj)) + 
+ (1 + cos^)^ sin {ojrf + e) t— 

-{l-cos^^)sm{{uJrf + e)t + 2ip)+ (13) 

+ (1 — cos^)^ sin [urf — e) t— 

- (1 - cos^ ^) sin {{urf -e)t + 2^)] . 
Fig. 4 shows the Fourier spectra of signals given by Eq. [T3] under the same conditions as in 

Figs. 2 and 3. For the random rf phase, the absorption signal has three comparable oscillating 

components with frequencies Ui and Ui ± \/bj\^ A'^^g (Fig. 4c). For the rf phase i/j = 0, the 

sidebands are smaller than those at the random rf phase by the factor A^s/ a/w2 + ^'bs (Fig- 

4a). When we use = 7r/2, the component with frequency ui vanishes and the sidebands are 

comparable to those at the random rf phase (Fig. 4b). Note that the high-frequency sideband is 

always more intensive than the low-frequency one. 




Figure 4: Fourier spectra of the absorption signals in 
the singly rotating frame shown in Fig. 2 by a dot 
line. 

Upon the rotating wave approximation {Abs — ?■ 0), it follows from Eq. [13] that for ip = only 
the component with frequency Ui remains. At the same time, for both ip = n/2 and the random 
rf phase, the intensities of the sidebands are equal. The equalization of sidebands can be used 
to indicate the validity of the rotating wave approximation. On the contrary, their asymmetry 
reveals the effect of the counter-rotating component of the rf field. Such asymmetry was observed 
in the dressed Rabi oscillations using the EPR experiment with random rf phase [5]. 



Note that, upon the resonant monochromatic interaction, the Mollow triplet is formed by 
the transitions between the dressed states of the ground and excited levels of the initial qubit. 
Similarly, at the doubly resonant bichromatic interaction, the triplet in the singly rotating frame 
is formed by the transitions between the split states of both the ground and excited levels of the 
initial qubit (Fig. 1). 

We now provide the expression for the absorption signal in the frame described by the density 
matrix P5(t) = e~*^5*p5(0)e*^^*, where ^5(0) is given by Eq. [TTJ In this doubly rotating frame 
(DRF), the absorption signal can be written as follows: 

vr^{t) = ((l|p5(t)|2) - (2|p5(t)|l)) /2i = (cos^cos^/^sinet -sin^/'cosrf) /2. (14) 

According to Eq. [TU the absorption signal in the doubly rotating frame is caused by the 
transitions between spin states dressed simultaneously by the mw and rf fields. At the exact 
resonance {ui = uJrf), the signal for ip = is smaller than the signal for ip = it/ 2 by the factor 
Abs/ \/i^2 + ^B5- If ^BS 0, the signal for = disappears. In this case, for = 7r/2, the 
absorption signal oscillates with the Rabi frequency U2- So, for ^ = 0, the absorption signal 
is fully due to the counter-rotating component of the rf field and its amplitude is proportional to 
the value of the Bloch-Siegert shift A^^. 

In conclusion, we have studied the evolution of spin qubits at the double resonance (wq = 1^, 
Ui = Urf) with a bichromatic field, consisting of transverse (high-frequency) and longitudinal 
(low-frequency) components. We have found that the double "dressing" of an initial qubit by the 
bichromatic field forms four new qubits with a smaller energy splitting, giving rise to multiphoton 
transitions. In the laboratory frame, three triplets correspond to the transitions between states of 
these qubits. The transition amplitudes depend strongly on the phase of the low-frequency field. 
The counter-rotating component of the low-frequency field causes the asymmetry of sidebands 
in the triplets. After taking into account this component, the errors in operations with qubits on 



dressed states in the strong-field regime are minimized. The types of operational multiphoton 
transitions can be selected by the choice of the rotating frame: one triplet, (cji, u)i± e), can be 
observed in the singly rotating frame, and only the transition at the frequency e is realized in 
the doubly rotating frame. 
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